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ABSTRACT 

The consequences of structured flows continue to be a pressing topic in relating spectral 
data to physical processes occurring in massive star winds. In a preceding paper, our group 
reported on hydrodynamic simulations of hypersonic flow past a rigid spherical clump to explore 
the structure of bow shocks that can form around wind clumps. Here we report on profiles of 
emission lines that arise from such bow shock morphologies. To compute emission line profiles, 
we adopt a two component flow structure of wind and clumps using two "beta" velocity laws. 
While individual bow shocks tend to generate double horned emission line profiles, a group of 
bow shocks can lead to line profiles with a range of shapes with blueshifted peak emission that 
depends on the degree of X-ray photoabsorption by the interclump wind medium, the number of 
clump structures in the flow, and the radial distribution of the clumps. Using the two beta law 
prescription, the theoretical emission measure and temperature distribution throughout the wind 
can be derived. The emission measure tends to be power law, and the temperature distribution 
broad in terms of wind velocity. Although restricted to the case of adiabatic cooling, our models 
highlight the influence of bow shock effects for hot plasma temperature and emission measure 
distributions in stellar winds and their impact on X-ray line profile shapes. Previous models 
have focused on geometrical considerations of the clumps and their distribution in the wind. Our 
results represent the first time that the temperature distribution of wind clump structures are 
explicitly and self-consistently accounted in modeling X-ray line profile shapes for massive stars. 

Subject headings: Stars: early-type - Stars: massive - Stars: mass-loss - Stars: winds, outflows - 
X-rays: stars 
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1. Introduction 

The subject of X-ray production in massive star 
winds continues to be an evolving field of study. 
The superionization seen at UV wavelengths of 
OB stars were best explained by a model that 
had a source of X-rays in the winds (Cassinelli, 
Castor, & Lamers 1978; Cassinelli & Olson 1979). 
Initial observations by the Einstein observatory 
made the important discovery that essentially all 
O stars were X-ray sources (Harnden et al. 1979; 
Seward et al. 1979). A key finding to emerge from 
these early observations is that the observed X- 
ray luminosities are roughly correlated with the 
bolometric luminosities as Lx ~ 10~^Lbo1 (e.g., 
Cassinelli et al. 1981). Additional more extensive 
studies confirmed the relationship (e.g., Berghofer 
et al. 1997; Naze et al. 2011), although the basis 
of the relationship continues to be a point of in- 
vestigation (e.g., Owocki & Cohen, 1999; Owocki 
et al. 2011). In addition, the majority of OB stars 
display soft X-ray emissions with temperatures 
kT < 1 keV (e.g., Berghofer et al. 1996; Giidel 
& Naze 2009); 

Two pictures for the X-ray emission from hot 
stars arose: one involving a coronal zone at the 
base of a cool wind (Cassinelli & Olson 1979) and 
one involving shocks that form by line-driven wind 
instabilities (Lucy & White 1980; Lucy 1982). The 
coronal model as the sole source of the observed X- 
ray emission was quickly ruled out based on analy- 
ses of the earliest higher spectral resolution obser- 
vations using the Solid State Spectrometer (SSS) 
on the EinsteAn observatory. Cassinelli & Swank 
(1983) found that the predicted large X-ray opti- 
cal depths expected for a base coronal source of 
X-rays were incompatible with the observed SSS 
spectra. They further suggested that these winds 
consist of many shock fragments to explain the 
lack of significant X-ray variability. 

Studies of X-ray emissions from OB stars have 
focused primarily on exploring the wind driven in- 
stabilities (or line de-shadowing instability, LDI) 
as a process of producing a distribution of wind 
shocks (e.g., Owocki, Castor, & Rybicki 1988). A 
detailed picture of the expected X-ray production 
from these wind shocks was given by Feldmeier 
(1995), and Felmeier et al. (1997) showed that a 
wide range of temperatures could be produced in 
a planar shock front. 



We are now in an era of high spectral resolution 
X-ray astronomy with a few dozen massive stars 
having been studied in long pointed observations 
(e.g., Walborn, Nichols, & Waldron 2009). Bet- 
ter quality data have led to a host of new ques- 
tions concerning the physics of X-ray generation 
in massive star winds (e.g., Waldron & Cassinelli 
2007; hereafter WC07). Most of the X-ray line 
emission is clearly formed within the winds. A 
triad of lines from He-like ions (forbidden, inter- 
combination, and resonance or "fir" lines) provide 
direct information about the formation radius of 
X-ray line emission (Kahn et al. 2001; Waldron & 
Cassinelli 2001; Leutenegger et al. 2006). Super- 
giant winds typically show that the lower energy 
ion stages such as O vii tend to form near or above 
10 i?*; intermediate energy ions (e.g., Ne IX and 
Mg Xl) form deeper at « 3 to 8 i?*; and high en- 
ergy ions such as Si xiii and S xv form relatively 
close to the star (< 2i?*). Waldron & Cassinelli 
(2001) suggested that these differences in depths 
could perhaps be explained from considerations of 
wind absorption effects, since the cool wind opac- 
ity scales as K oc A^. Thus winds are more trans- 
parent at shorter wavelengths (higher energies). 
Waldron & Cassinelli (2001) also noticed that the 
location of line formation for the He-like ions ap- 
peared to correlate with the respective radii of op- 
tical depth unity for the X-ray photoabsorption 
(c.f., Cassinelh et al. 2001; Miller et al. 2002; Os- 
kinova et al. 2006; WC07). The conclusion is that 
hot plasma is spatially distributed in the wind 
flow. 

One surprising result from the high resolution 
X-ray spectroscopy data is the general symme- 
try of broad lines and the frequent absence of 
line profiles with significantly blue-shifted peak 
emissions. It had been expected that the wind 
X-ray lines would be generally broad yet skewed 
to the short wavelength (or "blueward") side of 
the line. This skewness is a consequence the fact 
that for an expanding wind, the column depth 
of photoabsorption to the flow on the far side of 
the star is larger than on the near side, result- 
ing in differential attenuation between the red and 
blueshifted hemispheres (MacFarlane et al. 1991; 
Ignace 2001; Owocki & Cohen 2001). So the ob- 
servation of frequently symmetric and unshifted 
lines was unexpected for the massive winds of 
OB supergiants. For these stars WC07 found that 
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« 60% of emission lines are broad with a mean 
line-width (HWHM) of 0.3 to 0.5 of the wind ter- 
minal speed (Voo), and in excess of 75% of the lines 
have line-shifts that lie within ±0.2Voo of line cen- 
ter. 

As suggested by Waldron & Cassinelli (2001), 
the simplest way to account for the rather sym- 
metric and unshifted X-ray line profiles is that 
the wind is more optically thin to X-rays than 
suggested by the mass-loss rates. A variety of 
models have emerged to explain the line symme- 
try problem by studying wind clumping and wind 
porosity effects. Clumping in dense Wolf-Rayet 
winds has been known for many years (Moffat 
et al. 1988), and there is direct evidence of clump- 
ing among some 0-stars (e.g., Lepine & Moffat 
2008). Clumping can be categorized as ranging 
from micro-clumping (e.g., Hillier 1991; Hamann 
& Koesterke 1998) to macro-clumping (e.g., Feld- 
meier et al. 2003; Brown et al. 2004; Oskinova 
et al. 2004, 2006, 2007; Owocki & Cohen 2006), 
or a mix of the two. Micro-clumping explicitly 
assumes all clumps are optically thin at all wave- 
lengths, which need not be the case for macro- 
clumping. 

Reductions in the mass-loss rate M by a fac- 
tor of 10 or more appeared to be supported by 
FUSE observations of Pv lines from several hot- 
stars (FuUerton et al. 2006). Although this would 
certainly make the winds more thin to X-rays, this 
severe reduction in M can be eliminated either by 
accounting for wind "macro-clumping" (Oskinova 
et al. 2007) or by including the effects of XUV ra- 
diation in reducing the fractional abundance of Pv 
(Waldron & Cassinelli 2010). 

Other models that have been proposed to ex- 
plain the symmetry of the X-ray lines include 
the effects of resonance line scattering on line 
shapes (Ignace & Gayley 2002), and there is sup- 
port in one case where such effects are applicable 
(Leutenegger et al. 2007); two-component wind 
structures where the polar wind component is im- 
peded by surface magnetic structures (MuUan & 
Waldron 2006); and models requiring magnetic 
fields and coUisionless shocks (Pollock 2007). 

Resolved X-ray lines have served as an impe- 
tus to understand more accurately the nature of 
the hot plasma component. And encoded within 
these detailed X-ray emission line shapes is the re- 
quired information both about the formation pro- 



cess of the line (i.e., the density and temperature 
which determines the emissivity) and the vector 
velocity field. Although these various approaches 
have certainly had successes in trying to decode 
these line profiles, there remain open questions 
about understanding the temperature and emis- 
sion measure distributions, and the radial loca- 
tion of hot plasma formation and maintenance. In 
particular, previous considerations of line profiles 
from clumpy/porous winds (e.g., Owocki & Cohen 
2006; Oskinova et al. 2007) have focused on issues 
of clump geometry (pertaining to photon escape) 
and clump distributions (pertaining volume filling 
factors), but these have not self-consistently in- 
cluded temperature distributions implied by the 
structures themselves, as for example in planar 
shocks. Even smooth wind considerations have 
been geometrical in nature (e.g., Owocki & Co- 
hen 2001; Ignace & Gayley 2002). The models 
presented here have the benefit of self-consistently 
including the detailed temperature distribution of 
the shocked structures, within the context of the 
assumed model. 

The underlying model for clump bow shock 
structure was presented in Cassinelli et al. (2008; 
hereafter Paper I), who considered the shape, tem- 
perature, and density of bow shocks that form 
around wind clumps. In this second paper of 
the series, we are explicitly interested in the line 
profiles that form from these bow shock struc- 
tures and how features of the clump bow shock 
paradigm may contribute to understanding the ob- 
served shapes of massive star X-ray line profiles. 
In section [2] the results of Paper I are briefly re- 
viewed. In section [3l line profiles are calculated 
and discussed for individual clump bow shocks, 
emphasizing the diversity of line shapes that can 
result. Section |4] describes line profiles that from 
an ensemble of clumps, including the limitng case 
of many clumps and the case of a discrete en- 
semble of randomly placed clumps. Section [5] 
presents concluding remarks about these results 
and needed future areas of stdy. The Appendix 
details considerations of the temperature distri- 
bution in the wind for our model prescription. 

2. Model Description 

Our model calculations of X-ray emission line 
profiles produced by a wind distribution of clumps 
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and their associated bow shock structures are 
based on simulations discussed in Paper I. Our 
results apply to the hypersonic limit, namely that 
the Mach number is high ( > 10), which is an 
excellent description of the situation in a mas- 
sive star wind where the terminal speed Voc ~ 
1000 km s^^ and the gas thermal speed is w 100 
times smaller. This means that the bow shock 
structure - its shape, and its relative density and 
temperature distributions arc largely indepen- 
dent of the actual Mach number (e.g., Hayes & 
Probstein 2004). 

However, since line profile calculations require 
detailed information on the actual velocity field 
of a large number of wind distributed clumps, we 
need to establish the distribution of the line-of- 
sight (LOS) velocities as seen by an observer. The 
simulations of Paper I were carried out in the rest 
frame of a rigid and spherically symmetric clump 
where a fast moving wind with plane parallel sym- 
metry and constant density sweeps across the face 
of a clump. The plane-parallel approximation ap- 
plies when the the radius of the clump, Rd, is 
small compared to the wind clump radial location, 
r. To compute synthetic line profile shapes the 
vector velocity field for the bow shock found in 
Paper I (see their Fig. 4), which is accomplished 
through relative velocity vector defined by 

AV(r) = [Vw(r)-V,i(r)]r (1) 

where Vw is now the preshock ambient (or inter- 
clump) stellar wind speed at the site of the clump, 
Vci is the clump speed, and both are measured rel- 
ative to the stellar rest frame. We are assuming 
that Vw and Vd are purely radial and thus func- 
tions only of r. Note that the magnitude AV is 
the same in both the stellar and clump rest frames. 

Our calculations apply to both cases of a clump 
moving faster or slower than the ambient medium 
because the simulations are conducted in the rest 
frame of the clump. As long as the relative veloc- 
ity (AV) between the clump and the surrounding 
gas medium is hypersonic, the same bow shock 
structure results. The only practical difference for 
line profiles is in which direction the bow shock 
opens with respect to the star center. If the clump 
is moving radially outward faster than the wind, 
the stagnation point will be ahead of the clump, 
and the bow shock opens toward the star. Similar 



scenarios have been discussed by Guo (2010) and 
Waldron & Cassinelli (2009). If the clump is slow, 
then the geometry flips by 180°, like in the appli- 
cation to T Sco for infalling clumps (Howk et al. 
2000). 

An evaluation of the line profile shape requires 
knowledge of the LOS Doppler shifts toward an 
observer in a specified direction. Evaluating the 
Doppler shifts requires the introduction of several 
coordinate systems. We assume that the stellar 
wind is structured but spherical in a time aver- 
aged sense. With no special direction, the stellar 
and observer coordinate systems are chosen to be 
conincident. Cartesian coordinates {X, Y, Z) are 
introduced, with associated standard spherical co- 
ordinates (r, 1?, (p), where the polar angle is mea- 
sured from the +Z axis. The observer is located 
along the +Z axis. 

For each clump we adopt the cylindrical coor- 
dinate system (tu, (f>, z) used in Paper I, where the 
z-axis in cylindrical coordinates coincides with the 
symmetry axis of the bowshock, and w is the cylin- 
drical radius. The clump center corresponds to 
n7 = and z = 0. 

In addition, since we envision these clumps as 
moving radially from the star, the symmetry axis 
of the bow shock is also the radial line from the 
star center to the clump center. Hence, a key con- 
dition inherent in these two coordinate system def- 
initions is that z • f = ±1 is maintained for all 
wind distributed clumps, where the sign indicates 
whether the bow shock opens away from the star 
{+) or toward the star (— ). 

The specific points of our model are discussed 
in the following sections. We first start with a 
brief review of the bow shock properties found in 
Paper I as modified by using a relative velocity. 

2.1. Overview of the Bow Shock Proper- 
ties 

2.1.1. Geometry 

The numerical simulations of Paper I showed 
that the shape of the bow shock can be well- 
described by the form 

^— ^=a(— ) (2) 
with a = 0.35, m = 2.34, and zq = — l.lQiJcb 
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hence a shape not far from a parabola. The bow 
shock apex forms at a distance of 0.19i?ci above 
the clump surface. 

In addition to the bow shock shape, it was 
demonstrated that the derivative of the bow shock 
shape (i.e., the position-dependent tangent) is the 
key parameter in determining the velocity, tem- 
perature T, and emission measure EM distribu- 
tions along the bow shock surface. In Paper I we 
defined this derivative as 

-^ = 5(tn) = am(n7/i?ei)"-\ (3) 

It is convenient here to introduce an angle a that 
is related to the curvature of the bow shock, with 

tana = g(cc7). (4) 

Note that in Paper I, we had defined this angle as 
Ai but here prefer to use a. 

2.1.2. The Velocity Field 

One of the major findings from Paper I was 
that the EM was dominated by the immediate 
post-shock gas. So an "On The Shock" (OTSh) 
approximation was introduced, whereby the den- 
sity N ^ temperature, and velocity relevant to X- 
ray observables are described by conditions along 
the bow shock surface, and thus by the geometry 
described in the preceding section. Hence, the ve- 
locity field in the rest frame of the clump that is 
needed to synthesize line profiles is known analyt- 
ically at every point along the bow shock for the 
known surface geometry. 

The simulation was based on the assumption 
that the radius of the clump i?ci ^ r, so that 
the incident interclump wind flow was essentially 
plane parallel. Thus AV = AV z in the rest 
frame of a clump itself (c.f.. Fig. 4 of Paper I). 
We introduce the unit vectors n as the outward 
normal to the shock and I as a unit vector par- 
allel to the shock in the direction away from the 
apex. The jump conditions for a strong oblique 
shock were applied to derive the velocity com- 
ponents perpendicular and parallel to the shock 
front. Using primes to denote velocities in the 
clump rest frame, the post-shock velocity compo- 
nents are given by 



V'p,± = --|AF|cosan (5) 
1 1 , . , 

V'p.ii = |AV^| sinaf (6) 
= , ^ \AV\l 

where ISV is from equation ([T|). Note that /S.V can 
be positive or negative depending on whether the 
clump is traveling faster or slower than the inter- 
clump wind, but the perpendicular velocity com- 
ponent is always away from the bow shock sym- 
metry axis, and the parallel component is always 
downwind of the apex position. 

The postshock velocity fleld is a function of 
clump's radial distance from the star via IS.V 
and also location along the bow shock implicitly 
through the curvature factor g{vj). The total post- 
shock flow speec|3 anywhere along the bow shock 
is 

As expected, the postshock velocity takes on a 
value of Vp = Ay/4 for a head-on collision at 
the bow shock apex where (7 = 0. Far downstream 
of the apex, the speed approaches Vp — AV^ for 
g 1 as the shock becomes extremely oblique. 

In order to determine Doppler shifts of the X- 
ray emitting material, the postshock velocity fleld 
must be evaluated in the observer frame. To ac- 
complish this, it is convenient to express the vec- 
tor flow in the clump system as components of z 
and w which can readily be transformed to the 
star/observer system. This requires a standard 
rotation of coordinate systems from h and I to w 
and z, which is given by 

fi — w sin a =F £ cos a 
I = w cos a ± z sin a 

The distinction in sign for the z-component is re- 
lated to whether the clump moves faster or slower 

^We point out that eq.[7]corrects a typo appearing in eq. (10) 
of Paper I 
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than the interclump wind. If slower, then the bow 
shock opens away from the star and the upper sign 
is used; if faster, then the bow shock opens toward 
the star, and the lower sign is used. 

With these preceding expressions, the post- 
shock velocity in the rest frame of the star (un- 
primed system) becomes 



Vp 



[3 



1 



%m.2aw ± -(1 + 3sin^ a) z 
4 



+Vd z, 



(8) 
(9) 



where V4 > is always true, but > for the 
case that the clumps are slower than the inter- 
clump wind and < for clumps that are faster. 

The observed velocity shift for flow at a point in 
the wind \s vz — —Z - Vp, using lowercase "w" to 
signify that the velocity is for the observer. Car- 
rying out the dot product yields 



The path length downwind of the apex can be 
expressed as a function of impact parameter w for 
known bow shock geometry z{'uj). The postshock 
temperature, T, along the shock is found to be 



T{m) 



1 



(13) 



The power-law form of equation ([3]) provides gim), 
and the temperature distribution reduces to 



T 



= (1 + 9')-' = 



0.67 



Rci 



2.68 



(14) 

The line profile calculation also requires an 
emission measure distribution. A result of Pa- 
per I was that the differential emission measure 
{DEM = dEM/dT) is a monotonic power-law 
function of temperature and thus location along 
the bow shock with distance from the apex. One 
can conveniently parameterize the distribution 
with 



Vz = — \Vc\ + V:i\ cos -d + cos ip sin (10) 

Owing to axial symmetry of the bow shock, V4 
and 14 are functions only of distance from the bow 
shock apex along its surface. The velocity field 
reduces to Vp ~ 14 near the bow shock apex and 
Fp w 14 in the far outer wings of the bow shock. 

2.1.3. The Temperature and Emission Measure 
Distributions 

With the OTSh approximation, the runs of T 
and EM along the shock are monotonic functions 
of path length from the shock stagnation point. 
The peak temperature at the bow shock apex of 
an individual clump is Ta with a value given by 



Ta = ^'-^^V' 
16 k 



= 14 MK 



f— ) 
V0.62/ 



Ay 



1000 kms" 



(11) 
'(12) 



where in the latter expression we have evaluated 
the constants assuming a fully ionized gas with 
solar abundances. 



dEM _ EMo{r) 
~ TA{r) 



T 

tI 



-7/3 



(15) 



where the mapping between w and T is made 
through the factor g using equation The 
emission measure scaling parameter EMo is given 
by El 



SMo = 5.1x10^^ cm"^ 



(16) 

where we have assumed a strong shock such that 
the postshock number density is Np = ANw, 
where Nw is the preshock interclump wind num- 
ber density, and and have both been 
scaled to the values used in the simulation. From 
Paper I, we found that — 0.1i?ci well approxi- 
mates the simulation results. Our model does not 
predict the evolution of i?ci through the fiow. To 
compare most easily our results with those of pre- 
vious works, we adopt a scaling of R^^ oc r^. Im- 
plicit then is that the emission measure varies as 



^Eq. (22) in Paper I is missing a factor of (m — which 
leads to a slightly smaller scale factor for EMo as compared 
to the Paper I result. 
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EMo oc (rVw) ^ in form, like a smooth wind with 
a constant fihing factor of hot plasma (e.g., Ignace 
2001). 

2.1.4- The Interclump and Clump Velocity Dis- 
tributions 

The critical parameter of the bow shock model 
is the relative velocity defined in equation ([T]) since 
it determines the magnitude of both the postshock 
X-ray temperature structure and the velocity field. 
Our approach is to assume a two-component wind 
model where the ambient wind and clump entities 
follow different velocity laws. We adopt the com- 
monly used standard /3 velocity law prescription 
defined as 



V{r) = 



1 



(17) 



with Voo the terminal speed and 6 < i?, so that 
the radial wind speed is non-zero at the wind base 
taken to be the stellar radius i?*. The value of b 
sets the initial flow speed Vbj with 



1 



b 



/3 



(18) 



In all cases considered, the b parameter will be 
fixed at the same value in both the interclump and 
clumped wind velocity laws. In order to limit the 
number of free parameters for our exploratory in- 
vestigation, we also assume that both components 
achieve the same terminal speed. 

The smooth interclump wind component will be 
described by the parameter jSw- Throughout this 
paper, we adopt = 1 as typical of OB star wind 
solutions. Then j3c\ represents possible clump ve- 
locity distributions. Note that Pc\ > Pw implies 
slow moving clumps relative to the interclump 
flow; < Pw corresponds to fast moving clumps. 

2.2. Scalings from the /3-Law Prescription 

The distributions of Ta and the DEM can 
be derived from the two /? approach. We intro- 
duce the convenient velocity normalization w = 
Vw/Voo- With the terminal speeds for the clump 
and interclump flows the same, the relation be- 
tween the two components' velocity laws are w = 
w^j"', where Wc\ — K^i/Kxj- Thus the velocity jump 
of equation ([T]) becomes 



(19) 



This relation can be used to find Ta [r) , which pro- 
ceeds as follows. 

Equation (|12p along with the preceding expres- 
sion gives 



Ta ^ 



^ ^"^g AT/2 

16 k 



3 t^TTiH ,,2 



16 



w 



w 



1 - w^^' 

(/3ci-l) 



-1) 



(20) 

}(?1) 

(22) 



where Tiim is implicitly defined as the highest pos- 
sible temperature in our model of outfiow that oc- 
curs for a velocity jump that is equal to the wind 
terminal speed. 

Now the maximum hot plasma temperature 
Tmax in the wind model can be determined. In the 
velocity coordinate of the interclump fiow, T^ax is 
achieved at a critical value Wc as given by 



l/(/3ci-l) 



which in the clump velocity becomes 



Wcl. 



Pel 



/3ci/(/3ci-l) 



(23) 



(24) 



The radial location of T,„ax is at a corresponding 
critical radius value of Tc, with 

rc = — (25) 

1 - Wc 

The value of T^ax is determined by just two pa- 
rameters: the value of /3ci and the wind terminal 
speed via Tum, as given by 



T„,ax = Tun (/3cl - 1)' 



(26) 



Figure [T] shows the distribution of Ta in terms 
of the maximum possible temperature Tii„i with 
different curves for different values of /3ci. This 
is plotted against the normalized velocity of the 
interclump wind in the upper panel, and against 
the normalized velocity of the clumps in the lower 
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panel. The curves range from (3ci — 2 (lowest 
curve) to Pd = 8 (highest curve) in integer val- 
ues. As ^ci increases, Tmax shifts to progressively 
higher velocities of the interclump wind but lower 
velocities for the clump flow. Values of Ta at dif- 
ferent velocity locations are at the level of a tenth 
to a few tenths of Tum- For typical massive star 
wind speeds of 1000-3000 km s^^, Tn^ has values 
of 10-100 MK. 

3. Line Profiles for An Individual Clump 

Before developing emission line profiles for 
clumped winds, it is instructive first to con- 
sider the emission line shape arising from a single 
clump. As an example case, we consider a clump 
at a location of 2i?* that follows a /3ci = 3 velocity 
law. The velocity jump is AV « 0.4^00- Figure [2] 
demonstrates the diversity in profile shapes for 
this single clump when it is located at different 
positions around the star, as given by the angle 
illustrated by the inset. The abscissa is the LOS 
observer velocity shift wz = vz/Voo- Note that 
the profiles have been normalized to have unit 
area. Values of = 0°, 30°, 60°, 90°, 120°, 150°, 
and 180° were considered as labeled. For this 
figure both stellar occultation and absorption of 
X-rays by the clump itself are ignored, and the in- 
terclump wind is taken to be completely optically 
thin to X-rays. 

Except for = 0° and 180° which are clumps 
that lie along the LOS to the star, the profiles tend 
to be double-peaked and asymmetric. One excep- 
tion is when a clump is at = 90°; the profile 
is still double peaked but also symmetric since it 
lies in the plane of the sky with the star center. 
Generally, the double-horn shape is a consequence 
of the complex velocity field in the bow shock. 
The shapes become more nearly single-peaked as 
they approach the LOS to the star center. This 
is because the observer views the bow shock ex- 
actly along its symmetry axis. In conclusion, for a 
clump at 1? = 0° and 180°, only the component 
contributes to observed Doppler shifts, but for a 
clump located at z9 = 90°, only the Vt^ component 
contributes. 

For the profiles of Figure [U emission from the 
bow shock contributes from the peak tempera- 
ture Ta down to an imposed minimum of 0.5 MK, 
which we use as a low temperature cut-off for hot 



plasma X-ray production. However, real lines form 
only over a restricted temperature range with con- 
sequences for the line shape. Consider a hypo- 
thetical line that forms between 2 and 3 MK. For 
a bow shock with Ta = 10 MK, this line would 
arise spatially from an annular band centered on 
the symmetry axis of the bow shock and offset 
from its apex. Consequently, realistic lines that 
form over different temperature ranges will tend 
to have different shapes, because they sample dif- 
ferent portions of the postshock velocity field. 

Figure[3]illustrates this effect through the use of 
simple temperature cut-offs. The different curves 
are for line emission with different low temperature 
thresholds Tiq. Below Tio the emissivity is zero; 
above it the emissivity is independent of T. In 
this example clumps are placed at = 90°. The 
profiles becomes progressively broader as the lower 
temperature cut-off increases, with values of Tio = 
0.1 MK (blue), 0.3 MK (green), 1.0 MK (red), and 
3.0 MK (black). 

To understand the growing line width with in- 
creasing Tio, recall that the EM of a clump is dom- 
inated by the low temperature gas. For a clump 
a.t d — 90°, the bow shock is viewed perpendic- 
ular to its symmetry axis. Only V,^ components 
of the postshock velocity field contribute to ob- 
served Doppler shifts. With the lowest tempera- 
ture gas found furthest downwind of the bow apex, 
where the velocity vector is more nearly tangent to 
our LOS, Vcj tends to be relatively small. Lower 
speed flow in the bow shock is to be found closer 
to the apex; however, this flow has a relatively 
larger component in the m direction because of 
the greater curvature, with higher LOS Doppler 
shifts resulting at the bowhead. But the bow- 
head is exactly where the hottest plasma is to be 
found. Thus raising Tio means that the bowhead 
region increasingly dominates the line formation, 
typically leading to a broader line for the given 
geometry. 

Two final comments. First since increasing Tjo 
restricts the contributing volume, higher values of 
Tio also lead to weaker lines for a given clump. 
This is not apparent from Figure [3] because each 
profile is normalized to unit area. Second, the 
stagnation point at the bowhead is the hottest 
gas and has intrinsically very low speed flow in 
the clump rest frame. If Tio were to approach the 
value of Ta the profiles would actually narrow, a 
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limit not reached in the exampes of Figure |3l 

4. Line Shapes from an Ensemble of 
Clumps 

Although it is important to understand the 
emission profile from an individual clump bow 
shock, stellar winds are understood to be highly 
structured from many lines of evidence (e.g., Lupie 
& Nordsieck 1987; Hillier 1991; Mofatt & Robert 
1994; Lepine k Moffat 1999, 2008; Oskinova, Feld- 
meier, & Hamann 2004; Owocki & Cohen 2006; 
Prinja & Massa 2010; Muijres et al. 2011). Within 
our framework, this means there is more than one 
clump. Foremost is the basic observation that X- 
ray emissions from single massive stars are not 
highly variable. Although there is suggestive evi- 
dence of line variability (e.g., Nichols et al. 2011; 
Hole & Ignace 2012), in terms of bandpass lumi- 
nosities, OB stars are typically variable at the level 
of 10% or less (Cassinelli & Swank 1983; Berghoe- 
fer & Schmitt 1994; Berghoefer et al. 1996, 1997). 

Since we know that the observed X-ray emis- 
sion from these stars arises from a wind distribu- 
tion of X-ray sources, we now need to consider 
an ensemble of clump bow shocks for producing 
synthetic line profiles. We recognize the intrin- 
sic time-dependent nature of the problem which, 
in principle, requires a full radiation hydrodynam- 
ics approach (e.g., Dessart & Owocki 2003, 2005). 
Since a goal of this paper is to present an ini- 
tial analysis of line shapes arising from bow shock 
structures, such a detailed approach beyond the 
scope of this paper. Our basic premise is that 
observed emission lines reflect a time averaged 
wind flow. High energy resolution X-ray spec- 
troscopy of high signal-to-noise requires exposure 
times ranging from 50 to 200 ks. By contrast, the 
characteristic flow time in a massive star wind is 
R*/vac ~ 1 — 10 ks. This means that a typical mas- 
sive star X-ray spectrum is formed over multiple 
flow times, which tends to average over structural 
variations that are stochastic. 

4.1. The Limit of Many Clumps 

Having considered emission profiles from an in- 
dividual clump in section [3l here we consider the 
opposite extreme of many clumps, which we refer 
to as the effectively "smooth" limit. It is imagined 
that large numbers of clumps are uniformly dis- 



tributed in radius and orientation about the star 
to achieve strict spherical symmetry. Certainly, 
approximate spherical symmetry is consistent with 
low limits on the net continuum polarizations in 
O stars (McDavid 2000; Clarke et al. 2002). Po- 
larization of unresolved sources is related to devi- 
ations of a circumstellar envelope from spherical 
(e.g.. Brown & McLean 1977; Brown, Ignace, & 
Cassinelh 2000). 

The idealized smooth limit has value in estab- 
lishing a reference baseline of models that can be 
used for interpreting line profiles from winds with 
different degrees of structuring and time varying 
effects. (The latter will be treated in a separate 
paper.) Additionaly, the smooth limit allows for a 
derivation of the DEM from the wind as a whole. 
We begin this process by referring to equation ()15|) 
that describes the DEM of a single clump. The pa- 
rameters EMo and Ta are themselves functions of 
clump location via equations and (ITB)) . Both 
of these are functions of radius (or equivalently ve- 
locity) alone. Consequently, every clump in a shell 
of radius r will have exactly the same DEM. The 
global DEM for the wind will consist of integrating 
contributions provided by every shell. 

Working in the velocity coordinate of the in- 
terclump wind Vw, the total wind DEM is given 

by 



fdEM\ 



^^C^) dEM dN 



Vi(T) 



dT dVw 



dVw, (27) 



where dN/dV represents the clump distribution 
in terms of radial velocitjU. In the limiting case 
of many clumps uniformly distributed in velocity, 
dN/dV is a constant. To find the total DEM, 
the integration proceeds only over shells where the 
temperature T is high enough to produce X-ray 
emission. 

As described previously, there is a maximum 
temperature located in the wind at radius with 
corresponding normalized velocity Wc- The de- 
scription of apex temperatures is thus double val- 
ued with radius. Plus, any given shell will have 



•^For example, Sundqvist, Owocki, & Puis (2011) show a 
clump filling factor as a function of radius, both as inferred 
from observations and deduced from model simulations. 
Their single peaked curve from H* to large r corresponds 
to a bell-shaped distribution in dN/dV from Vq to Voo- 



9 



a range of temperatures from Ta down to a lower 
cut-off value. Clearly, a particular value of T will 
only be found in a shell if Ta{Vw) > T. It is 
this condition that is used in equation (P7| for the 
limits of the integrand. 

Using results from the preceding section and 
equations (fT2|) and (fT6|) . equation ([27l) becomes 



f dEM \ 
\ dT I 



oc 





— 7/3 / 


J lim / 




{w^ 





1 — Ul 



(28) 



where wi and ?i;2 represent the velocity interval 
between which T is achieved, and F is a correction 
factor for stellar occultation. The latter is given 

by 



2 




(29) 



Note that the integral is over a fixed T, hence the 
temperature dependence T~'^/^ can be factored 
out of the integral. Thus the integral that remains 
represents a temperature-dependent modification 
to the power law for a single clump. The inte- 
grand is not overly complex, but the integration 
limits wi{T) and W2{T) tend not to be analytic. 
(See App. B for solutions of wi and W2 in the spe- 
cial cases of (id = 1/2, 2, and 3.) 

Figure |4] displays the results of calculations for 
the total DEM of the wind at even values of /3ci 
from 2 to 12. The DEM is plotted logarithmically 
against T/Tmax- The curves have been shifted to 
a zero value at the lowest temperature used. The 
results all lie very close to each other with depar- 
tures from a —7/3 slope occuring only at higher 
values of T, as emphasized by the dotted line for a 
power-law of —7/3 slope. Different /3ci values yield 
the —7/3 slope at low T because in our model the 
cooler X-ray emitting plasma is to be found essen- 
tially throughout the wind. The slight steepening 
toward larger T becomes a downturn as T^ax is 
approached, because only the hottest components 
are severely restricted in radial locale. 

The range of line profiles that result in the 
smooth limiting case are displayed in Figures [5] 
and m A stellar wind terminal speed of Voo = 
2500 km s~^ is adopted, as before. For all cases 



there is a minimum radius (or velocity) in which 
hot X-ray emitting plasma is to be found. For the 
upper panels, the radius is r,„in = l.li?*, and for 
the lower panels, it is r-cnin — 1.5i?,. The two left 
panels are profiles that result for fast clumps with 
Pel = 0.5; the two right panels are for slow clumps 
with /3ci — 3. 

The line emission is assumed to be optically 
thin (hence no resonance scattering effects). The 
model line luminosity as a function of relative ve- 
locity shift wz — vz/Voo is calculated by 



-Xwz) 



A{T)^e-^dT (30) 



where A(r) is the temperature dependent line 
cooling function, r is the photoabsorption optical 
(see below), and the integral is carried out over 
the unocculted volume. For initial calculations we 
assume simply that A(T) is a constant. We also ig- 
nore the variation of ion fraction with T, implicit 
in the DEM factor. In effect, these illustrative 
model line profiles are meant to sample the full 
DEM distribution. The inclusion of T-dependence 
for line cooling function and the effects of ioniza- 
tion balance should result in a more diverse set of 
line profiles. 

The different colored curves in Figures [S] and 
ini are for different levels of wind photoabsorption. 
The assumption is that the clumps are small com- 
pared to other scales in the problem so that the 
photoabsorption optical depth r is approximated 
from a LOS integration through a smooth inter- 
clump wind. As such, the photoabsorption optical 
depth to the clump position is the same for all 
points along the bow shock. 

The optical depth is calculated following Ignace 
(2001), with 



r(r,i?) 



sin-i? 



1 



[{bsm^')/{rsm-&)y 

(31) 

where is the optical depth scale to the base of 
the wind at i?, . Generally, this scale is related to 
the wind mass-loss rate, abundances, and energy 
of the particular line transition in question. 

In Figures [5] and |6] profiles for values of r, = 
0,0.5,1,2,4, and 8 are calculated. The effect of 
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increasing is to make the profiles increasingly 
asymmetric with emission peaks of progressively 
higher blueshifts. The black profile is for a line 
with T* = 0; magneta corresponds to the case of 
T* = 8. Note that with no photoabsorption, the 
black curve displays a central "fiat-top" indicative 
of Vniin = V'(''min), modulo the cffcct of Stellar 
occultation. 

For the range of photoabsorptive optical depths 
used, the blueshifted peaks all lie below about half 
of terminal speed. The line widths actually de- 
crease slightly for low r*, but then increase with 
larger values of r*. Ultimately, at large optical 
depths, the blueshift of the peak emission and the 
line width are not sensitive to the value of rmin. 

4.2. The Case of Discrete Clumps 

Relaxing the assumption of a smooth distri- 
bution of clumps, a discretely structured fiow is 
now considered. Models are based on a random 
number generator to place clumps throughout the 
wind from which the X-ray emission line profiles 
are computed. Let s be a random number in the 
range of to 1. Individual clumps are sprinkled 
in a uniformly random way about the star. The 
ith clump will have angular coordinates given by 
ipi = 2iTs and /Xj = —1 -|- 2s, where = cosi^i, 
and of course two distinct random numbers s are 
used to set the two seperate coordinates for a given 
clump. 

Radial placement requires a different approach. 
In our model clumps exist exterior to the photo- 
spheric level, i?*, out to infinite distance, in princi- 
ple. The ID radiative hydrodynamic simulations 
indicate that the formation of strong shocks oc- 
curs primarily at low and intermediate radii in the 
flow (e.g., Feldmeier, Puis, & Pauldrach 1997), ba- 
sically where the velocity gradient is reasonably 
strong. Structure can persist and evolve out to 
fairly large radius. As a way of capturing the flavor 
of this scenario, we choose to space clumps such 
that they are statistically uniformly distributed in 
radial velocity Vw- In relation to the preceding 
section, this means that dN/dV now becomes a 
uniform probability distribution to be sampled in 
the range of Vq to Voo- 

The relationship between a random number s 
and the corresponding velocity for that value is 
given by: 



Naturally, this distribution is highly non-uniform 
in radius. On average half of the clumps lie at 
Vw > 0.5 Voo, and the other half lie below that 
speed. For (3^ = 1, this means that half the 
clumps lie beyond r — 2i?*, and half lie inte- 
rior. One can easily incorporate different distri- 
butions dN/dV, either as an exploration of pa- 
rameter space or to match known clumping prop- 
erties of a particular source. The manner in which 
X-ray producing structures are placed in velocity 
space influences the line profile shape. Our choice 
of dN/dV as uniform is merely a convenience for 
purposes of illustration. 

Figure [7] shows examples of emission lines for 
different clump ensembles. The wind photoab- 
sorption optical depth is set to a low value of 
T* = 0.1. All profiles have been normalized to 
unit area and so no vertical scale of flux is pro- 
vided. The six panels labelled (a)-(f) correspond 
to different numbers of clumps A/'ci with 4 in (a), 
8 in (b), 16 in (c), 32 in (d), 64 in (e), and 128 
in (f). The black curves are the intrinsic profiles 
of the model calculation, whereas the red curves 
are convolved by a Gaussian to simulate the ef- 
fect of instrumental smearing from finite spectral 
resolution. 

It is important to note that the number of 
clumps contributing to a given profile is generally 
less than the value of A/'ci. This occurs for a couple 
of reasons. First, we adopt a threshold tempera- 
ture of 0.5 MK for gas to contribute to the line. If 
the apex value Ta is less than the threshold, then 
all the gas in the bow shock of that clump is also 
less than the threshold. The threshold eliminates 
those clumps that are very near the photosphere 
and very far away, where AV is too small to gen- 
erate the requisite temperatures for X-ray emis- 
sion. The second reason is that some clumps are 
occulted. 

With A/'ci on the order of several tens and 
higher, the convolved profiles are reasonably sym- 
metric (but not exactly so). Of course the extent 
of blueshifted peak emission and line width is a 
function of photoabsorption optical depth. 

We have not properly dealt with the fact that 
there is generally a broad range of temperatures 
across the bow shock. The emission lines of Fig- 
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ure [7] still adopt a temperature independent line 
emissivity as was used for the effectively smooth 
wind case of section 14.11 A temperature depen- 
dent emissivity should be included when fitting 
observed line profiles for specific sources. 

5. Summary and Conclusions 

Paper I of this series presented results of a 
hydrodynamic simulation for purely adiabatic 
cooling with a plane-parallel hypersonic fiow im- 
pinging upon a rigid spherical obstacle in the 
rest frame of that obstacle. The simulation was 
conducted under the assumption that individual 
clump structures are much smaller than the radius 
at which they are located. In that paper the fiow 
and temperature structure were described, and 
two quite interesting simplifications were empha- 
sized. First, it was found that the DEM followed 
a power-law form. Second, the emission measure 
was to be found primarily in a thin "sheath" of 
postshock volume. Thus Cassinelli et al. (2008) 
introduced the on-the-shock approximation, or 
"OTSh" , whereby the bow shock geometry deter- 
mines the T and DEM distributions necessary for 
computing observables. 

In this second paper, we adopt the OTSh to 
model X-ray emission lines that would arise from 
an individual bow shock and from an ensemble 
of bow shocks. This follows on a long string of 
papers to explain the unexpected observed X-ray 
line profile shapes from a number of massive stars 
in terms of structured fiows, based on fragments 
of planar shocks (Oskinova et al. 2004) or porosity 
arguments (Owocki & Cohen 2006). 

An individual clump tends to produce an asym- 
metric double- horned emission profile that is offset 
from line center, depending on its radial and lat- 
eral location around the star from the observer. 
Evidence indicates that massive star winds are 
characterized by large numbers of clump struc- 
tures. To model the line shapes from an en- 
semble of clumps, we adopted a parametric two- 
component flow approach using two wind /3-laws: 
one for the interclump wind flow and one for the 
clump fiow. The distinction in /3-laws leads to 
radius-dependent velocity jumps that govern the 
temperature range of the bow shocks. Of particu- 
lar interest is that this approach yields a number 
of semi-analytic relationships for the T and DEM 



distributions throughout the fiow, which in prin- 
ciple are properties that can be tested against ob- 
servations (e.g., Waldron & Cassinelli 2009; Guo 
2010); 

Using this construction, emission line profiles 
were calculated in the "smooth" limit of many 
uniformly distributed clumps and for the case of 
a discretely structured fiow. As expected, peak 
emission of the lines are a function of the de- 
gree of photoabsorption. The bow shock paradigm 
yields line shapes that are somewhat symmetric at 
modest photoabsorption optical depths of a few, 
where the infiuence of rmin on the line shape can 
no longer be perceived. In contrast to a uniform 
distribution of clumps, the discrete case leads to 
profiles with spikey features; however, these are 
much too narrow to actually resolve with current 
instrumentation. Using a simple temperature cut- 
off approach, we also find that profile widths can 
depend on the temperature interval of line forma- 
tion. 

All of these results represent a promising start- 
ing point for tailored analyses of individual ob- 
jects, for calculating spectral energy distributions, 
and for investigating X-ray variability. Previ- 
ous efforts have focused primarily on geometri- 
cal considerations for explaining X-ray line pro- 
files shapes observed from OB stars, in the form of 
discrete clumps, clump distributions, and/or fill- 
ing factor considerations. Our results explicitly 
include temperature distributions throughout the 
wind fiow, which is a forward step in X-ray line 
profile synthesis modeling. 

In closing we remind the reader that our ap- 
proach has relied on simulations that adopt purely 
adiabatic cooling for the bow shocks. We have 
begun new simulations of clump bow shocks that 
include radiative cooling. The advantage of adia- 
batic cooling and hypersonic fiow is that the fiow 
geometry is independent of the Mach number. 
In situations where radiative cooling is needed, 
the results will depend on the density and the 
apex temperature achieved. Consequently, the 
bow shock structure will no longer have a "uni- 
versal" form; thus, greater complexity is the cost 
of greater realism. Preliminary results with ra- 
diative cooling suggest that the power-law DEM 
in temperature derived in Paper I persists at the 
hottest temperatures, but shows a flattening to- 
ward cooler temperature gas where radiative cool- 
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ing dominates. In the futiire we will include the 
results of these new simulations along with realis- 
tic temperature-dependent line emissivities to fit 
the line profiles of high resolution X-ray lines from 
massive star winds and to study time variable ef- 
fects of X-ray emissions. 
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A. Appendix: Temperature Intervals for /3c\ — 1/2, 2, and 3 

To calculate the total DEM from a uniform distribution of many clumps in a wind, it is necessary to find 
the integration limits wi(T) and W2(T) in eq. (1291) . Based on the preceding section, this amounts to a root 
finding exercise involving the following relation (see eq. |22j): 

w^-' -w + Vi = 0, (Al) 

where t = T/Tn^ and l3w = 1 is assumed. The function is double-valued for all /3ci 7^ 1- Note that the 
clump /3ci can be larger or smaller than the interclump value. Here solutions are given for three cases where 
the roots are analytic or semi-analytic. 

A.l. Case of /Jci = 2 

The equation to be solved is 



w"^ - w + Vi = 0. (A2) 



The roots have with values of 



wi,2 = 1^1 \/l - 4Vi. (A3) 
The maximum temperature occurs at imax — 1/16 for which wi ~ W2 = 0.5. 

A.2. Case of ;9ci = 1/2 

The equation to be solved is 

Vt = Q. (A4) 
With the change of variable x"^ = w, the condition can be recast as 



2 x + Vt^O, (A5) 



which is the same quadratic expression for /3ci — 2. The roots wi,2 for the case /3ci ==1/2 are simply the 
square roots of the solutions from the = 2 case. The maximum temperature still occurs at tmax = 1/16, 
which in velocity is now wi = W2 = 0.25. 

A.3. Case of /3ci = 3 

The expression to be solved is 



Vi = 0. (A6) 



This cubic has three real roots; however, one of those is negative and not physical. There are standard forms 
for the roots; here we use the trigonometric version. An angle 7 is introduced as defined by 



cos 7 = - VV^max, ( A7) 

4/27. Then the roots become 

W2 = ^ cos(7/3), (A8) 



15 



and 



Wi 



1- 



1 

77 W2- 



(A9) 
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Fig. 1. — A plot of bow shock apex temperatures Ta for clumps located at different positions in the wind. 
The temperature is normalized to Tum (see text). The upper panel shows location in terms of the interclump 
wind velocity; lower is for the clump velocity. Curves are for different ^ci values, ranging from 2 (lowest 
curve) to 8 (highest curve) in integer intervals. 
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Fig. 2. — The inset (top center) shows the location of a clump at angle i) around the star from the observer's 
axis. The plot shows example emission line profiles, all normalized to have unit area, for individual clumps 
located at the indicated orientations. In each case the clump is at the same radius, and so all profiles have 
the same apex temperature Ta- Solid curves are for clumps on the nearside of the star; dashed are for ones 
on the far side. 
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Fig. 3. — Similar to Fig. [5] but now profiles are for clumps only at i9 = 90° and with different temperature 
intervals. The emissivity is taken to be constant within the temperature range of Tio up to Ta, with 
Tio = 0.1, 0.3, 1.0, and 3.0 MK from the most narrow line (blue) to the broadest one (black), respectively. 
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Fig. 4. — A logarithmic plot of the intrinsic total DEM against temperature in the smooth wind limit. 
Temperature is normalized to Tmax- Curves are for /3ci values of even integers between 2 to 12, inclusive. 
The dotted line is for a —7/3 power law as would apply to a single clump. The curves have been shifted to 
have the same value at the lowest temperature for ease of comparison. Despite the wide range of /3ci values, 
similar overall DEM distributions result. 
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Fig. 5. — Illustrative emission line profiles for the smooth limiting case. Hot plasma is assumed not to 
exist interior to r^-m, with a value of l.li?, for the upper panel and 1.5i?* for the lower one. These models 
are for fast clumps with /3ci = 0.5. Different colored curves are for different levels of interclump wind 
photoabsorption, with values of t* = 0, 0.5, 1, 2, 4, and 8 that lead to increasingly blueshifted lines. 
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Fig. 7. — Line profile shapes for an ensemble of clumps with — 0.1. Panels are distinguished by the 
number of clumps Afd used in the model, with (a) 4, (b) 8, (c) 16, (d) 32, (e) 64, and (f) 128 clumps. Model 
line profiles are shown in black; overplotted are red curves that include the effects of instrumental smearing 
are included. Finite spectral resolution is approximated by convolving model lines with a Gaussian that has 
cr = O.OSKo. 



23 



